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Abstract 

We give a method which generates sufficient conditions for instability of equilibria for circulatory 
and gyroscopic conservative systems. The method is based on the Gramians of a set of vectors whose 
coordinates are powers of the roots of the characteristic polynomial for the studied systems. New 
instability results are obtained for general circulatory and gyroscopic conservative systems. We 
also apply this method for studying the instability of motion for a charged particle in a stationary 
electromagnetic field. 
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1 Introduction 

Many physical phenomenons are modeled by second order Euler-Lagrange equations 

IdL dL 

dtdc? dq* V " [ ' 

where L is the Lagrangian function and Qi are the components of generalized forces. The linearized 
equation at an equilibrium point has the form 

Mq + A 2 q + A 3 q = 0, (1.2) 

where M is a constant symmetric matrix and A2, A3 are constant matrices. The above system appears, 
for example, in vibration theory and it is known as lumped-mass system (without external forces) (see 
PQ), where q is an n vector of time- varying elements representing the displacements of the masses. 
Assuming that M is invertible we multiply the equation ( |1.2[ ) by the matrix M _1 . Using the notations 
D and K for the symmetric parts, respectively G and C for the skew-symmetric parts of M _1 A2 and 
M~ l A% we obtain the following normal system, see pQ, [2] 

q+(D + G)q+(K + C)q = 0. (1.3) 

In literature, D is called the viscous damping matrix, G is the gyroscopic matrix, K is the stiffness 
matrix and C is the circulatory matrix. 

According to [T] we have the following classification of the normal systems of the form (1.3 1: 

(i) conservative systems when D = G = C = and K is positive definite; 

(ii) gyroscopic conservative systems or undamped gyroscopic systems when D = C = 0; 



* Corresponding author; West University of Timi§oara, Department of Mathematics, Bd. Vasile Parvan, No. 4, 300223 
Timispara, Romania; tel.: +40 726 424147, fax: +40 256 592316, e-mail: birtea@math.uvt.ro 



(iii) damped non-gyroscopic systems or passive systems when G = C = and D, K are positive definite; 

(iv) circulatory systems when D = G = 0; 

(v) systems with constraint damping when G = 0. 

In Section 2 we give a list of sufficient conditions for the existence of a complex root with non-zero 
imaginary part for a polynomial with real coefficients. We obtains these conditions using the Gramian of 
a set of vectors whose coordinates are powers of the roots for the given polynomial. We explicitly write 
three of these conditions, which turn out to be conditions involving only the power sums of the roots 
for the given polynomial. In the case of a characteristic polynomial of a matrix, the above mentioned 
conditions are written only in terms of norms and traces of the symmetric and skew-symmetric parts of 
the matrix. 

In Section 3 we give sufficient conditions for the instability of circulatory systems. These conditions 
involve only the norms and traces of the stiffness matrix and the circulatory matrix. We recover 
a sufficient condition previously found by Bulatovic in [3]. We also give an example of a circulatory 
system, which shows that the three sufficient conditions presented in this section give different instability 
regions. 

In Section 4 we give sufficient conditions for the instability of gyroscopic conservative systems in- 
volving norms and traces of the stiffness matrix and the gyroscopic matrix. We apply these conditions 
for the case of a charged particle acted by the Lorentz force. We also describe the regions of instability 
for this system. 



2 Sufficient conditions for complex roots with non-zero imagi- 
nary part 

Let Q(a) = a n + aia™^ 1 + ■ ■ ■ + a n be a n-degree polynomial with real coefficients. We denote by 
{«!, . . . , a n } the set of roots of the polynomial Q. We search for sufficient conditions in order to have 
at least one complex root with non-zero imaginary part. 
Consider the vectors: 

v = (!,..., 1); 
vi = (ai, . . . , a n ); 



v n - 1 = (a%- 1 ,...,o%T 1 ). 
We have the well-known result: if all roots of the polynomial Q are real, then for any < i\ < ■ ■ ■ < 



ik < n — 1 the following inequality holds 
Gram^, ...,v ih ) := 



> 0. 



< v ik ,v n > .. 

Lemma 2.1. If there exist < i\ < ■ ■ ■ < ik < n — 1 such that 

Gram(u il , . . . ,v ik ) < 0, 
then the polynomial Q has at least one complex root with non-zero imaginary part. 

n 

The power sums of the roots are real numbers denoted by Sk ■— a ii where k 6 N. We notice that 

i=l 



Gram(w , ■ ■ -,v n -i) = 



n si 
si s 2 

Sn—l $n 



Sn-1 



Sln-1 



■2 



Applying Lemma 2.1 for the simple particular cases Gram(i>o, vi), and Gram(uo,i>2) an d Gram(i>i, V2) 
we obtain the following sufficient conditions. 



Proposition 2.2. // any of the inequalities hold 
(i) ns 2 < s\; 
(ii) ns4 < s|; 
(iii) S2S4 < s§, 

then the polynomial Q has at least one complex root with non-zero imaginary part. 

Considering Gramians of higher order, one can obtain similar sufficient conditions. We notice that 
the inequalities (i) and (ii) hold when s 2 < 0, respectively S4 < 0. Also, inequality (iii) holds, in 
particular, when S2 < and S4 > or S2 > and S4 < 0. If P(A) = Q(X 2 ), then Proposition 



2.2 



gives 

sufficient conditions for the existence of a complex root with strictly positive real part of the polynomial 
P. 

All the power sums Sk can be expressed recurrently, using Newton formulas, in terms of the coeffi- 
cients of Q only (see [1]): 



Sk 
Sk 



- aiSfc-i 
- ais fe _i 



+ ak-isi + kak = 0, 1 < k < n 

+ dk-iSk-n+i + a n s k - n = 0, k > n. 



The above recurrences and Lemma |2.1| give sufficient conditions for the existence of complex roots 
with non-zero imaginary part exclusively in terms of the coefficients of the polynomial Q. More precisely, 
we have the following result. 

Proposition 2.3. The polynomial Q has at least one complex root with non-zero imaginary part if: 
(i) 



(ii) 



n = 2 
n — 3 
n > 4 



n>2: n(a 1 — 202) < a 1 ; 

a\(a\ - 4a 2 ) < 0; 

a\ + 6ai<23 + a\ < 4a^a2; 

n(a 1 — Aa\a2 + 4aia3 + 2a\ — 4a4) < {a\ — 2a 2 ) 2 ] 



(iii) 



n = 2 
n — 3 
n > 4 



a\{a\ — 4a 2 ) < 0; 



a\a\ + 100,1020,3 < 2a^fa3 + ia 2 



a\a\ + \Qa\a2az 



9a 



80204 < 2a^ja3 



3; 
4a^a4 



4a^ + 9ai 



For the case when Q is the characteristic polynomial of a matrix M £ M„> 



i.e. 



Q(a) 



det(a/„ — M), conditions (i), and (ii) and (iii) of Proposition 2.2 can be interpreted only in terms of the 
symmetric and skew-symmetric parts of the matrix M. Any matrix M can be uniquely decomposed as 
a sum of a symmetric part and a skew-symmetric part, M = M s + M a , where M s = h(M + M T ) and 
M a = l(M-M T ). 

In what follows, we will use the well-known equalities: Tr(A + B) = Tv{A) + Tr(B), and Tr(AB) = 
Ty(BA) for any A, B £ 3vl nxn (R.) and Tr(AB) — for any symmetric matrix A and any skew-symmetric 



3 



matrix B. Furthermore, we have ||A|| 2 = Tr(AA T ) for any A = [ay] € M„ X „(R), where ||A|| 2 := a %- 

id 

Consequently, ||M S || 2 = Tr(A/ 2 ) and ||AT Q || 2 = — Tr(M 2 ). By a direct computation we have 
Sl = Tr(M) = Tr(M s ); 

s 2 = Tr(M 2 ) = Tr(Af 2 + M 2 a + M s M a + M a M s ) 

= Tr(M 2 + M 2 ) = \\M s \\ 2 -\\M a \\ 2 ; 
s 3 = Tr(Af 3 ) = Tr((Af s + M a f) 

= Tr(Af 3 + Af s Af 2 + M 2 s M a + M s M a M s + M a M 2 s + M a M s M a + M 3 + M 2 M S ) 

= Tr(Af 3 ) + 3Tr(Af s Af 2 ); 
s 4 = Tr(Af 4 ) = Tr((M s + M a f) 

= Tr(Af 4 ) + Tr(Af 4 ) + 4Tr(M 2 M 2 ) + 2 Tr((Af s A/ a ) 2 ) 

= ||Af 2 || 2 + ||A/ 2 || 2 - 4\\M s M a \\ 2 + 2 Tr((Af s Af Q ) 2 ). 

The next theorem is a completion of a result given in [3] and it follows by substituting the above 
expressions for 81,82,83, S4 in Proposition |2.2| 

Theorem 2.4. If one of the inequalities hold 

(i) n(||A/ s || 2 -||Af a || 2 ) <Tr 2 (M s ); 

(to) n(||Af 2 || 2 + ||M 2 || 2 -4||Ar s Af Q || 2 +2Tr((A/ s A/ a ) 2 )) < (||M S || 2 - ||Af a || 2 ) 2 ; 

(Hi) (\\M S \\ 2 -\\M a \\ 2 )(\\M 2 \\ 2 + \\M 2 \\ 2 -4\\M s M a \\ 2 + 2Tr((M s M a ) 2 )) < (Tr(Af 3 ) + 3 Tr(A/ s A/ 2 )) 2 , 
then the matrix M has at least one complex eigenvalue with non-zero imaginary part. 

For n = 2, if Tr(Af ) =^ 0, then the inequalities (i) and (ii) from the above theorem are equivalent. 
As we will show later, this is not the case for n > 3. 



3 Instability for circulatory systems 

First we will give sufficient conditions of instability for a circulatory system subject to potential (con- 
servative) and non- conservative positional (circulatory) forces 

q + /vq + Cq = 0, (3.1) 

where K is symmetric and C is skew-symmetric. In [3] it has been showed that if the circulatory forces 
are bigger, in a certain sense, than the conservative forces, then the equilibrium point q = 0,q = is 



unstable. The characteristic polynomial of the system (3.1) is given by P(X) = det(X 2 I n + K + C) . This 
polynomial contains only even powers of A. Denoting A 2 = a we obtain the characteristic polynomial 
Q(ct) = det(aJ„ + K + C) of the matrix M = —{K + C). The existence of a complex root a with 
non-zero imaginary part implies the existence of a complex root A for the polynomial P with strictly 



positive real part. Thus, instability of the system (3.1 1 is implied by the existence of a complex root for 



the polynomial Q(a). Applying Theorem 2.4 we obtain a different proof of a previous instability result 



presented in [3.: and two new instability criteria. 
Theorem 3.1. If one of the following inequalities hold 

(i) r0/;n(||^|| 2 -||C|| 2 )<Tr 2 (^); 

(ii) n{\\K 2 \\ 2 + ||C 2 || 2 - ±\\KC\\ 2 + 2Tt((KC) 2 )) < (\\K\\ 2 - \\C\\ 2 ) 2 ; 

(m) {\\K\\ 2 - ||C|| 2 )(||^ 2 || 2 + ||C 2 || 2 - 4||^C|| 2 + 2Tr((^C) 2 )) < (Tr(X 3 ) + 3Tr(XC 2 )) 2 , 
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then the equilibrium point q = 0, q = of the system (3.1 ) is unstable 



Remark 3.1. The inequality ns4 < s\ is equivalent with n ( a i) 2 < a f) 2 ■ Noticing that af, i = 

i=l i=l 

1,2,..., 7i, are the roots for the characteristic polynomial Q(a) — det(aJ n — M 2 ) we can apply the 
inequality (i) of Proposition 2.2 to the polynomial Q(a) = det(ce/„ — M 2 ). We obtain 



\\Ml\\ 2 >\\M 2 s \t 



1 



Tr 2 (M 2 ). 



An elementary computation gives us that M 2 = KG + CK and M 2 = K 2 + C 2 . Consequently, the 



inequality (ii) of Theorem 3.1 is equivalent with 

\\CK + KC\\ 2 > \\K 2 + C 2 \\ 2 



n 



\K\ 



\\C\\ 2 f 



A 



The inequalities S2 < or S4 < obviously imply the existence of at least one complex eigenvalue 
with non-zero imaginary part. As a consequence, we obtain the following sufficient conditions for 
instability: 

Corollary 3.2. If one of the following inequalities hold 
ft) 11*11 <» 

(ii) \\K 2 \\ 2 + \\C 2 \\ 2 + 2Tr((KC) 2 ) < A\\KC\\ 2 , 



then the equilibrium point q = 0,q = of the system (3.1) is unstable 



Other sufficient instability conditions for circulatory systems can also be found in [5], [B]. 

For n = 2, if Tr(K) 7^ 0, then the inequalities (i) and (ii) from the above theorem are equivalent. 
For n > 3 the conditions (i), and (ii) and (iii) give different instability regions, as the following example 
shows. Let 

K = 

k 

Condition (i) is equivalent with the inequality 3c 2 — 3fc 2 — 1 > 0, condition (ii) is equivalent with the 
inequality — (c 2 — fc 2 ) 2 + 14c 2 — 2fc 2 — 1 > and condition (iii) is equivalent with 4c 6 — 4fc 6 — 12c 2 fc 2 (c 2 — 
k 2 ) + 8c 4 — ik A + 4c 2 fc 2 + 4c 2 > 0. The Figure 1 below shows the three instability regions corresponding 
to conditions (i), (ii), respectively (iii). 



" 1 





" 







c 


" 





1 


k 


; c = 


— c 











k 


















Fi gure ll Instability regions using Theorem |3.l| the horizontal grid designates the instability region obtained with 
condition (i); the right-inclined grid designates the instability region obtained with condition (ii); the left-inclined grid 
designates the instability region obtained with condition (iii). 

The above example shows that conditions (ii) and (iii) provide two supplementary instability regions 
that cannot be determined by using condition (i). 
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4 Instability for gyroscopic conservative systems 



Next, we will give sufficient conditions of instability for a mechanical system subject to gyroscopic forces 
and potential forces 

q+Gq + ^q = 0, (4.1) 

where G is skew-symmetric and describes the gyroscopic forces and K is symmetric and describes the 
potential forces. Various results concerning the stability problem for gyroscopic conservative systems 
can be found in [7J, 0, [S], [10]. 



The matrix associated with the linear system (4.1 ) is given by 

On In 



A = 



-K 



-G 



and its characteristic equation is -P(A) := det(A 2 /„ + AG + K) = 0. Using the skew-symmetry of G and 
the symmetry of K we obtain that if A is a root of the characteristic equation, then —A is also a root of 
the characteristic equation. It follows that the above characteristic polynomial contains only powers of 
A 2 . After the substitution A 2 — a we obtained the reduced polynomial Q(a). We will apply Proposition 
2.2 and Proposition 2.3 for the polynomial Q(a) in order to obtain sufficient conditions for the existence 
of complex roots with non-zero imaginary part. Thus, we get sufficient conditions for instability of the 



equilibrium point q = 0, q = of the linear system (4.1 ) 



We denote by sj? the power sums of the roots of the polynomial Q(a). We also denote by the 



/i.S'o 



= 2 S 2 

a Q\2 



2s? 

Consequently, the condition (i) of Proposition 2.2 applied to the polynomial Q(a), 



power sums of the roots for the characteristic polynomial -P(A). The f ollowing equalities hold s 2 
and 



< (sj ) is equivalent with 2ns% < (s% ) , which is also the condition (ii) of the same proposition 
applied to the polynomial P (A). The inequality ns® < (sf) 2 implies that the polynomial Q(a) has at 
least one complex root with non-zero imaginary part and thus inequality 2ns ^ < (sf) 2 implies that the 
polynomial P(X) has at least one root with strictly positive real part, which guarantees the instability 

Tr(A 4 ). We have the following computations 



of (4.1). 



As before, = Tr(A 2 ) and s% 



A S = -{A + A 1 ) 



On ' 
Wn-K) () 



Wn-K) 



A n 



1 



(A - A T ) 



On 
\(In+K) 



\{In+K) 

-G 



We denote by B\ := \{I n — K), Pi '■— k{Jn + K). Consequently, we have 



At 
A 2 A 2 



B\ On 
On B\ 

B\ On 
On Bf 

-BM 



Al 
At 



-B% 



B 2 G 



GB 2 G 2 
— B 2 G 2 B'2 + B 2 



B\G 2 - B\B\ 



(A s A a ) 



(B1B2) 



{B x Bz 



Using the equalities Tr(G 2 B|) = Tr(S|G 2 ) = Tr(GB|G) = Tv(B 2 G 2 B 2 ) we obtain 

1 



Tr(G 4 ) + i (n + 4Tr(K) + 6Tr(K 2 ) + 4Tr(K 3 ) +Tr(K i )) 
8 

[Tr(G 2 ) + 2Ty{G 2 K) + Tr(G 2 K 2 )}; 



G 4 - G 2 B 2 - B\G 2 - GB 2 G + B\ 



Tx{K 2 



4Tr(A 2 A 2 ) = -^ 
2Tr((AsA a ) 2 ) = \ (n - 2Tr(X 2 ) + Tr(^ 4 )) 



Tr(G 2 ) - 2Tr(i^G 2 ) + Tr(K 2 G 2 ) - - Tr(K 4 ); 
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1 1 

- n - 
2 V 



- (n + 2Tr(K) + Tr(K 2 )) 



Using the equality Tr(G 2 K) = Tr(KG 2 ) we have 
Tr(A 2 s )+Tr(Al) 

2 Tr(K) + Tr(K 2 )) + Tr(G 2 ) 
-2Tr(^)-||G|| 2 ; 

Tr(At) + Tr(Ai) + 4Tr(^4 2 ) + 2 Tv((A s A a ) 2 ) 
= 2 Tr(K 2 ) + Tr(G 4 ) - 4 Tr(GifG) 
= 2 Tr(K 2 ) + Tr(G 4 ) + 4 Tr(G T KG) 
= 2\\K\\ 2 + \\G 2 \\ 2 + ATr(G T KG). 
Summarizing, we obtain the following instability result that was also achieved in with a different 



4 



argument. 

Theorem 4.1. If the following inequality holds 

2n{2\\K\\ 2 + \\G 2 \\ 2 + 4Ti-(G T KG)) < (2Tr{K) 
then the equilibrium point q ; 



|G|| 2 ) 2 



0, q = of the system (4.1 ) is unstable 



For n — 3 we will study the instability for the classical example of a charged particle in a stationary 
electromagnetic field. The particle is acted by the Lorentz force = <?E(x) + qx x B(x), where q 
is the electric charge of the particle, E(x) is the electric field, B(x) is the magnetic field and x is the 
position of the particle. From the stationarity assumption and Maxwell-Faraday equation we have that 
E(x) = V0(x), where <j)(x.) is the electric potential. The equation of motion is x = ^V^(x) + ^xxB(x), 
where m is the mass of the particle. At an equilibrium point (x e , 0) we have that V</>(x e ) = and the 
linearized equation at this equilibrium point is 



y = — Hess 0(x e )y H B(x e )y, 



(4.2) 



where B(x e ) is the constant skew-symmetric 3x3 matrix associated to the constant vector B(x e ). By 

a convenient choice of coordinates the matrix — B(x e ) can be rendered in the canonical form 

m 



in 



B(x e ) 



-c 
c 




We consider the potential function 0(x) = -x 

-— Hess(/>(x e ), the equation (4.2) is in the form (4.1), where 



6X2X3. Making the notations G 



-B(x e ) and 



K 



' k 





" 







c 


" 








k 


; G = 


— c 











A: 

















K 



with k = -52 and c, k ^ 0. The characteristic polynomial is P(X) = A 6 + (c 2 + fc)A 4 - fc 2 A 2 - k 3 . After 
substitution A 2 = a, the reduced polynomial becomes Q(a) = a 3 + (c 2 + k)a 2 — k 2 a — k 3 . Condition 
i) of Proposition 2.3 applied for the polynomial Q(a) is equivalent with the condition from Theorem 
4.1 and gives the inequality — (c 2 + k) 2 — 3k 2 > 0. This inequality does not give an instability region. 



Condition (ii) of Proposition 2.3 applied for the polynomial Q(a ) is equivalent with the inequality 
c 6 + 4c 4 fc + 10c 2 fc 2 + 6fc 3 < and condition (iii) of Proposition 2.3 applied for the polynomial Q(a) is 
equivalent with the inequalities k < and 2c 6 + 7c 4 k + 18c 2 k 2 + 8k 3 > 0. The Figure 2 below shows 



the two instability regions corresponding to conditions (ii) and (iii) of Proposition 2.3 applied to the 
polynomial Q(a). 

For n = 2 the characteristic polynomial can be expressed in terms of traces and determinants of the 



matrices which appear in the normal form (1.3). The study of the instability region (when the matrices 
D and G are also present) has been completely solved in [T5], [13], using a result from [T4"] . 
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Figure 2: Instability regions using Proposition |2. 3[ the right-inclined grid designates the instability region obtained with 
condition (ii); the left-inclined grid designates the instability region obtained with condition (iii). 

5 Conclusions 

In this paper we have used the Gramian technique for obtaining sufficient conditions for the existence 
of at least one complex root with non-zero imaginary part for a polynomial. Applying these results 
we have obtained sufficient instability conditions for circulatory systems and gyroscopic conservative 
systems. In this way, we recover the main theorems in [3] and Using other Gramians of second or 
higher order, one can obtain similar sufficient instability conditions. Gramians have also been used in 
studying the geometry of a higher order dissipation [15] . 
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